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The Bose-Einstein condensate (BEC) of a dilute gas of bosons is well described by the 
three-dimensional Gross-Pitaevskii equation (3D GPE), that is a nonlinear Schrodinger equa- 
tion. By imposing a transverse confinement the BEC can travel only in the cylindrical axial 
direction. We show that in this case the BEC with attractive interaction admits a 3D bright 
soliton solution which generalizes the text-book one, that is valid in the one-dimensional 
limit (ID GPE). Contrary to the ID case, the 3D bright soliton exists only below a critical 
number of Bosons that depends on the extent of confinement. Finally, we find that the 3D 
bright soliton collapses if its density excedes a critical value. Our results are obtained by 
using a nonpolynomial Schrodinger equation (NPSE), an effective one-dimensional equation 
derived from the 3D GPE. 

§1. Introduction 

In a dilute gas of Bosons at zero temperature practically all particles are in the 
same single-particle state of the density matrix p(r, r'; t).El)'cP This macroscopically 
occupied single-particle state is called Bose-Einstein condensate (BEC) and it is 
well described by a complex classical field ip{r,t) (order parameter or macroscopic 
wavefunction) whose Lagrangian density is given by 
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where U{r) is the external potential, g = AttH ag/m is the interatomic strength 
with Ug the s-wave scattering length, and the complex field is normalized to A^. By 
imposing the least action principle one obtains the following Euler-Lagrange equation 
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which is a nonlinear Schrodinger equalion (NLSE) and it is called three-dimensional 
Gross-Pitaevskii equationJap GPE).0)'i) 

In recent experiments eP'Q^ Bose-Einstein condensates have been trapped in quasi- 
ID cylindrical optical traps to study bright soliton solutions, which are shape- 
invariant matter waves with attractive interatomic interaction (a^ < 0). In this paper 
we derive analytical results for the BEC bright solitons by using a one-dimensional 
nonpolynomial Schrodinger equation (NPSE) we have recently obtained from the 3D 
GPE.i'B 
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§2. From 3D GPE to NPSE 

The cylindrical optical confinement of Ref. 3 and Ref. 4 can be described by a 
harmonic potential 

U{r) = ^mojl{x^ + y^) , (2-1) 

in the transverse direction of the BEC. This external potential suggests to map the 
3D GPE into an effective ID equation, which simplifies greatly the solution of the 
3D GPE. Our approach is to choose for the complex field the following variational 
ansatz 

i:{r, t) = N^'^<t){x, y, t- a{z, t)) f{z, t) , (2-2) 
where both (f) and / are normalized to one and (p is represented by a Gaussian: 

g 2CT(z,t)2 

y, t; a{z, t)) = . . (2-3) 
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Moreover we assume that the transverse wavefunction (j) is slowly varying along the 
axial direction with respect to the transverse direction, i.e. V'^c/) ~ ^^(p where 
= d'^ jdx^ + d'^ jdy^ . By using the trial wave- function in the Lagrangian density 
and after spatial integration along x and y variables the Lagrangian density becomes 
a function of o", / and /*. The Euler-Lagrange equation with respect a is given by 



=ai\/l + 2a.Ar|/|2 , (2-4) 



where a_L = yJTi/ [mu±) is the oscillator length in the transverse direction. Inserting 
this result in the Euler-Lagrange equation with respect to /* one finally obtains 
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This equation is a time-dependent non-polynomial Schrodinger equation (NPSE). 
Note that under the condition a^A^I/p <^ 1 one has cr^ = aj_ and NPSE reduces to 
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where the additive constant huj± has been omitted because it does not affect the dy- 
namics. This equation, called one-dimensional Gross-Pitaevskii equation (ID GPE), 
is the familiar nonlinear cubic Schrodinger equation. The nonlinear coefficient g' of 
this ID GPE can be thus obtained from the nonlinear coefficient g of the 3D GPE 
by setting g' = g/{2TTa\). Note that the limit asN\f\'^ <C 1 is precisely the condition 
for the one-dimensional regime where the healing length ^ = (STrasyo)"^/^ of the BEC 
is larger than aj_. In Ref. 6 and Ref. 7 we have verified that NPSE is very accurate 
in the description of cigar-shaped Bose condensates both in the 3D regime and in 
the ID regime. NPSE has been recently applied in the study of an atom laser. Ef 
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§3. 3D BEC Bright Solitons 



Under transverse confinement and negative scattering length (a^ < 0) a bright 
sohtonB' sets up when the negative inter-atomic energy of the BEC compensates the 
positive kinetic energy such that the BEC is self-trapped in the axial direction. The 
shape of this 3D condensate bright soliton can be deduced from NPSE. Scaling z in 
units of aj_ and t in units of uj^ , with the standard position 



f{z,t) = 0{z - ^t)e™(^-^*)e^(^'/2-A.)t ^ 
from NPSE one obtains the Newtonian second-order differential equation 
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where C, = z — vt and 7 = \as\N/a^. A simple constant of motion of this equation is 
given by 
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from which one can write 
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By imposing the boundary condition (p ^ for (" — > 00, which implies that E = 0, 
by quadratures one obtains the solitary bright-soliton solution written in implicit 
form 
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Moreover, by imposing the normalization condition to the function one has 

(1-^)^/^-1(1-^)^/^ + ^7 = 0. (3-6) 



The normalization relates the chemical potential /i to the coupling constant 7, while 
the velocity v of the bright soliton remains arbitrary. Note that in the ID limit 
(7(^2 <^ iliQ normalization condition gives /i = 1 — 72/2 and the bright-soliton 
solution reads 
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The above solution is the text-book ID bright soliton of the ID nonlinear (cubic) 
Schrodinger equation (ID GPE). 

From Eq. (3-6) it easy to show that for 7 > 2/3 there are no solitary- wave 
solutions. This is a remarkable result because, contrary to the 3D bright soliton the 
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widely studied ID bright soliton exists (and it is stable) at any 7. Thus BEC bright 
solitons exist only below a critical number 

N. = (3-8) 
3 \as\ 

of Bosons. Finally, we observe that Eq. (3-5) is well defined only for 7^^^ < 1/2; at 
7^2 = 1^2 the transverse size a of the BEC soliton is zero. Because is measured 
in units of l/aj_, it follows that it exists a critical axial density 

^ (3-9) 
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above which there is the collapse of the BEC bright soliton. This result is confirmed 
by the numerical integration of 3D GPE. 



§4. Conclusions 



It is well known that in one dimension soliton states exist for Bose condensed 
atoms with negative scattering length described by the one-dimensional Gross- 
Pitaevskii equation, that is the one-dimensional nonlinear cubic Schrodinger equa- 
tion. Such states, called bright solitons, are characterized by a time invariant shape. 
We have studied the existence of similar states in higher dimensions. In two and three 
dimensions in a flat potential stable bright solitons do not exist. cP However if an ex- 
ternal harmonic potential constrains the motion in two dimensions we have found 
that such 3D soliton states exist but only below a critical number Nc of bosons. The 
value of Nc depends on the transverse conflnement, more precisely it is proportional 
to the ratio of the transverse harmonic length to the scattering length. Moreover, 
we have found that 3D bright solitons collapse if their axial density is larger than 
l/(2as): for this value the transverse width of the bright soliton shrinks to zero. 
Finally, it is important to stress that our solitary-wave solutions can be called soli- 
tonic solutions because they have particle-like nroperties, e.g. collisional stability, as 
recently verified studying scattering processes. EP 
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